Vortex Multiplication in Applied Flow: a Precursor to Superfluid Turbulence 
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A surface-mediated process is identified in 3 He-B which generates vortices at roughly constant 
rate. It precedes a faster form of turbulence where inter- vortex interactions dominate. This precursor 
becomes observable when vortex loops are introduced in low-velocity rotating flow at sufficiently 
low mutual friction dissipation at temperatures below 0.5 T c . Our measurements indicate that the 
formation of new loops is associated with a single vortex interacting in the applied flow with the 
sample boundary. Numerical calculations show that the single- vortex instability arises when a helical 
Kelvin wave expands from a reconnection kink at the wall and then intersects again with the wall. 

PACS numbers: 67.40.Vs, 47.37.+q, 98.80.Cq 



In superfluid 3 He-B a hydrodynamic transition takes 
place below 0.6 T c from regular (laminar) vortex flow at 
high temperatures to turbulent flow at low temperatures 
1]. The transition occurs as a function of damping in 
vortex motion, the mutual friction dissipation, which in- 
creases roughly exponentially with temperature. If 
a bundle of vortex loops is injected in applied flow, be- 
low 0.6 T c they interact generating new loops in a rapidly 
growing tangle of vortices. The onset temperature of this 
turbulence depends on the number of injected loops. At 
low enough temperature even a single injected vortex ring 
leads to turbulence. This is surprising since turbulence is 
thought to result from the collective interaction of many 
loops. Also in numerical calculations one ring does not 
lead to turbulence when placed in uniform applied bulk 
flow. What is the explanation of this conflict? 

Liquid helium flow is generally contained inside solid 
walls and thus the interaction of the expanding vortex 
ring with the sample boundary has to be examined. Our 
measurements indicate that initially the injected ring 
generates new vortices while interacting with the con- 
tainer wall. This is observed as slow vortex formation 
which precedes the later more rapid turbulence and has 
a lower onset temperature. It supplies new vortices so 
that ultimately rapid turbulence will switch on at some 
location where the loop density has grown sufficiently. 

The existence of a slow precursor, which later esca- 
lates to rapid turbulence, is duplicated in numerical sim- 
ulations. They show that a helical Kelvin wave 0, Q 
expands on a single vortex when it becomes aligned suf- 
ficiently parallel to the flow. A growing Kelvin wave may 
then reconnect at the boundary, creating one new vortex 
as well as sharp kinks. These kinks excite new Kelvin 
waves [f| , starting a self-repeating process of vortex mul- 
tiplication. We call multiplication a process in which new 
vortices are formed as a result of the dynamic evolution 
of existing vortices. 



Experiment: — We measure with NMR the evolution in 
the number of vortices when vortex loops are introduced 
with different techniques in applied flow. Our sample of 
3 He-B is contained in a quartz tube of length d — 110 mm 
and inner radius R — 3 mm. The flow is created by ro- 
tating the sample around its axis with constant angular 
velocity fi. Initially the sample is vortex- free and the ap- 
plied flow arises from the difference of the normal and su- 
perfluid velocities, the counterflow velocity v ci = v n — v s . 
If viewed from the laboratory, the normal component is 
in solid body rotation, v n = ft x r, while the superfluid 
component is stationary, v s = 0. The maximum veloc- 
ity Vd m — QR is at the cylindrical boundary. Such a 
state is possible if i> c f m is maintained below a container- 
dependent critical value @- This requirement is here 
observed so that the growth in vortex number N(t) mon- 
itors vortex multiplication after injection. 

In the over-damped regime of vortex motion T > 
0.6 T c , the injected loops expand to rectilinear lines, con- 
serving their number. In their lowest energy state they 
form a central cluster where the N straight lines are 
packed with an areal density n = 2Q/k. Here k = hjlm^ 
is the circulation quantum. Outside the cluster the coun- 
terflow increases from zero to i> c fm = [1 — A '/ ' (irR 2 n)]£lR 
at the cylindrical boundary. Well within the under- 
damped regime T < 0.6 T c , the injection ultimately 
always leads to rapid turbulence. Its signature is to 
create close to the equilibrium number of lines, N < 
ir R 2 (2fl I k) . N can be deduced at the top and bottom of 
the long sample by measuring the NMR spectra with two 
independent spectrometers |l( . The technique is based on 
the strong dependence of the spectrum on v c f . A calibra- 
tion can be constructed experimentally or by calculating 
the order parameter texture and its NMR spectrum nu- 
merically 0]- 

In the onset regime T < 0.6 T c , the evolution after in- 
jection depends on the injection method. The highest 
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FIG. 1: The probability that vortex multiplication will start 
from a remnant vortex, when rotation is increased from f2j 
to plotted vs. temperature (bottom axis) and mutual 
friction parameter q — a/(l — a') (top axis). No vortices are 
formed above 0.5 T c . Below 0.45 T c vortex formation always 
leads to turbulence. At intermediate temperatures only part 
of the runs result in vortex formation. 



onset for rapid turbulence is measured by making use of 
the properties of the AB interface between a short sec- 
tion of magnetic-field stabilized 3 He-A and the remaining 
3 He-B in the long sample Q . The AB interface undergoes 
an instability when Q is increased across a well-defined 
critical value il c AB = 1.2 - 1.6 rad/s and a bundle of 
~ 10 closely spaced loops is tossed on the B-phase side. 
Measured in this way, the onset is at 0.59 T c (at 29.0 bar 
pressure). Its half width of 0.03 T c we interpret to reflect 
the variation in the number of injected loops. 

The lowest onset is measured by exploiting vortex for- 
mation following absorption of a thermal neutron in 3 He- 
B lj. This reaction heats a blob of ~ 100 (im diame- 
ter to the normal state. While the blob rapidly cools 
back to the ambient temperature a vortex ring with a 
diameter similar to the blob size is created. The ring 
starts to inflate and evolve, if Q exceeds a critical value 
51 C n = 1.4 rad/s. At 0.59 T c a single ring does not lead to 
turbulence, but at 0.45 T c well above 80 % of the neutron 
absorption events at £1 = 1.6 rad/s develop to turbulence. 

This comparison shows that an additional mechanism 
(requiring a lower value of damping) is needed to start 
turbulence from a single injected ring than when a bun- 
dle of many loops is used. We assume that in the latter 
case rapid turbulence is switched on immediately at the 
injection site, while in the former case a precursor mech- 
anism is required. In both cases at £1 ~ 1.4 rad/s the flow 
velocity is so high that rapid turbulence follows within 
seconds and no time is left to monitor a precursor. 

These two injection techniques can be accurately con- 
trolled, but here the flow velocity cannot be reduced to 
arbitrary low values. For this two new methods were de- 
veloped. Conceptually simplest is the case of the remnant 



vortex at £1 = 8]. When rotation is stopped, vortices 
are rapidly pushed owing to their mutual repulsion to the 
boundaries for annihilation. Only the annihilation time 
of the last one or two vortices becomes long below 0.5 T c , 
owing to the small dissipation a(T) and the long sample 
length. A single straight vortex parallel to the cylinder 
axis at a distance b = SR from the center survives for a 
time [2ttR 2 /(an)] [—In 5 — |(1 — 5 2 )], where the prefac- 
tor equals lh at 0.4 T c . In practice the annihilation time 
can be longer, since the sample and rotation axes cannot 
be aligned perfectly parallel and the last vortex is not 
straight. If the time t w spent at £1 = is shorter, then 
the sample contains a remnant vortex of complex shape. 

In Fig. n we plot the probability of observing rapid 
turbulence after the rotation has been increased to fif = 
1.1 rad/s with a remnant vortex in the sample. The 
abrupt temperature dependence in this figure is charac- 
teristic of a transition as a function of rapidly varying 
mutual friction. No conventional explanation in terms 
of critical velocities or vortex mill behavior is capable 
of producing such a steep jump. It is centered at 0.47 T c 
while the corresponding value with injection from the AB 
interface is at 0.52 T c (at 10.2 bar). In both cases the half 
widths are 0.03 T c . We interpret that these two transi- 
tions do not overlap and that the case of the remnant 
vortex in Fig. ^ exemplifies the transition in the single- 
vortex regime where a precursor is required. 

To reveal the precursor, we reduce the applied flow ve- 
locity, ie. the value of £!f. The result is shown in Fig. |2 
where we plot the growth in the vortex number N(t) with 
time. The new feature is the slow approximately linear 
increase in N (solid line), before rapid turbulence sets in 
(dashed lines). At flf = 0.6 rad/s the slow increase lasts 
for more than 200 s, generating approximately 1 vortex 
in 5 s until about 20% of the maximum number of vor- 
tices have been created. At larger fif the slow increase is 
shorter in duration, eg. at ilf = 1 rad/s rapid turbulence 
starts already after 30 s. The slow vortex generation at 
roughly constant rate N we identify as the precursor. 

Two further observations about the precursor follow 
from Fig. [2 (1) Vortex formation proceeds independently 
in different parts of the sample. At fif — 0.6 rad/s it takes 
more than 300 s for a vortex created at one end of the 
sample to reach the other end 0. Still, vortex formation 
at the top and bottom is observed to proceed at roughly 
the same rate. Thus the precursor is not localized (as 
would be eg. a vortex mill). 

(2) In Fig. a second method has been used to in- 
troduce vortices in flow, to start from a more controlled 
initial vortex configuration. The equilibrium vortex state 
at a low rotation of £!i = 0.05 rad/s is employed as start- 
ing point. Here a few vortices in the outer peripheral 
ring next to the cylinder wall are not rectilinear but con- 
nect to the wall (since the sample and rotation axes can 
only be aligned to within about one degree, see inserts 
in Fig. [21 • To appreciate the influence of these curved 



300 



'e 200 

O 



1 



100 



P = 
T=0.45T r 




l.Orad/s 













\v n -i 


s 





in 1 r 
0.75 rad/s 



0.6 rad/s 



a 




o 0.6 rad/s, top 
♦ 0.6 rad/s, bot. 
A 0.75 rad/s, top 
▼ 0.75 rad/s, bot. 
□ 1.0 rad/s, top 
► 1.0 rad/s, bot. 



-100 



100 200 
Time t (s) 



300 



400 



FIG. 2: Vortex formation with a slow precursor which sud- 
denly develops to rapid turbulence. The number of vortices 
N(t) is recorded at the top and bottom of the sample. The 
solid line is the average slow rate N while the dashed lines de- 
note the rapid turbulence. Initially the sample is in the equi- 
librium vortex state at fii = 0.05 rad/s with N ~ 37 vortices 
and a few vortices connecting to the cylinder wall (left insert). 
Rotation is then increased to a new stable value f2f, which is 
reached at t = 0. Three runs with different Sl{ are shown. 
During the ramp to fif, N ~ const while the flow builds up 
and compresses the vortices in a central cluster (right insert). 
The precursor is attributed to an instability of the single vor- 
tices, which extend across the counterflow region and end at 
the cylindrical boundary. 



vortices, the experiment was repeated differently. 

A cluster with only rectilinear vortices, but with less 
than the equilibrium number, can be prepared at higher 
temperatures and can then be cooled to T < 0.6 T c . As 
long as this cluster is separated by a sufficiently wide 
counterflow annulus from the cylindrical boundary, fi can 
be increased or decreased without change in N at any 
temperature down to our lowest value of 0.35 T c . If is 
reduced too much, the cluster makes contact with the 
cylindrical boundary, some outermost vortices become 
curved, and during a subsequent increase of fi, while 
T < 0.5 T c , the behavior in Fig. [21 is reproduced. There- 
fore to observe vortex multiplication at f2f, we conclude 
that a curved vortex connecting to the cylindrical wall 
is required. The number of these initially curved vor- 
tices increases with J7i and thus also the rate N increases 
with fli. Nevertheless, the behavior is also observed when 
Oi = and one starts from a remnant vortex. 

At a low value of f2f a curved vortex which connects 
to the cylindrical boundary spends a long time expand- 
ing axially along the sample. The vortex segment ad- 
jacent to the wall moves in the counterflow created by 
the rotation and reorients partially along the flow owing 
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FIG. 3: Simulation of the measurements in Fig. [21 The total 
number of separate vortices N(t) vs. time is followed in a 
cylinder of radius 3 mm and length 10 mm. The calculation 
is started from a vortex ring of 2 mm radius in the azimuthal 
plane. The ring is unstable in azimuthal flow and generates 
via the Kelvin wave instability 3, 4] tens of vortices in a rapid 
burst, which gives the configuration shown in the insert. After 
this initial burst slow vortex formation at constant average 
rate N starts. Here each new vortex is produced from Kelvin- 
waves expanding on an isolated vortex (as shown in the insert) 
which is blown up to ring-like shape and then reconnects at 
the boundary. The later rapid growth in vortex number is 
dominated by inter- vortex interactions. 



to its self-induced velocity. It is thus expected to be- 
come unstable with respect to the formation of Kelvin 
waves 0, . The expanding waves may then reconnect 
with the wall and generate new separated loops. This 
interpretation for the linear precursor in Fig. [2] explains 
qualitatively its distinct features such as its abrupt tem- 
perature dependence or its low threshold velocity. To 
analyze the precursor in more detail we examine vortex 
dynamics numerically. 

Simulation: — The velocity «l of a vortex line element 
is determined from 3] 



v s + as x (v n - v B ) - a's x [s x (v n - v B )] . (1) 



Here s is a unit vector parallel to the vortex line element 
and a and a' are the mutual friction parameters mea- 
sured in Ref. 0]. From the form of this equation we can 
expect that the solutions for i>l can be classified by the 
parameter q = a/(l — a'). We calculate the evolution of 
the vortex configuration from Eq. JU, with v s obtained 
from the Biot-Savart law and the boundary conditions 
derived from an additional solution of the Laplace equa- 
tion 0, 0| . We apply background flow of v n in different 
geometries, including rotating flow in a cylinder and flow 
in a pipe with a uniform or parabolic velocity profile. 

We find that at low vortex density new vortices are 
generated only from expanding Kelvin waves which in- 
tersect with a wall. The growth of these waves depends 
strongly on the orientation of the vortex segment with 
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FIG. 4: As a model of the precursor, the collision of a vortex 
ring with a plane wall (at z — 0) is calculated. A ring of 
0.5 mm radius is initially above the wall, with the plane of 
the ring tilted by n/10 from the x = plane. In the frame of 
reference of this figure, uniform flow of the normal component 
at v n = 3mm/s is applied in the x direction. The different 
contours show the ring at 0.05 s intervals. The sharp kinks 
at the wall reconnections induce Kelvin waves on the loop. 
These are able to grow in the applied flow only on one of the 
two legs formed in the reconnection (here on the right). The 
largest wave reconnects with the wall and a new loop is then 
separated. The conditions for the growth of Kelvin waves are 
fulfilled in such collisions only in the under-damped temper- 
ature regime: The insert on the top shows the probability 
that a new vortex is created when the original ring is initially 
placed at z = 1 mm with random orientation. This probabil- 
ity depends on the ratio q = a/(l — a'), shown on the top 
axis, rather than on a or a' separately. 



nation of this difference is an interesting physical ques- 
tion. Self-sustained vortex multiplication has been pre- 
viously demonstrated in highly inhomogeneous bulk flow 
of the normal component [12j. When walls are present 
we find that the inhomogcncity of the applied flow is not 
essential to maintain continuous growth in N. Never- 
theless, the flow geometry affects the probability to start 
self-sustained vortex generation from a single seed loop, 
which is, for instance, larger in a circular pipe than in a 
rotating cylinder (at the same value of q). 

Conclusions: — Our results demonstrate that in the 
under-damped regime of vortex motion intrinsic vortex 
formation has to be described in an ideal superfluid as a 
sequence of multiple processes. It starts with the nucle- 
ation of the first vortex, is followed by surface-mediated 
multiplication of more vortex loops, and it finally goes 
over in rapid turbulence when inter-vortex interactions 
become possible. These individual steps have been diffi- 
cult to separate in earlier work. This is because in super- 
fluid 4 He with strong surface pinning the solid surfaces 
are covered with a plentiful source of remanent vortices 
which easily start rapid turbulence when flow is applied. 
A second reason is that experimentally oscillating flow 
has been easiest to achieve and there the accumulation 
of a vortex tangle occurs differently, as was recently re- 
ported from measurements with a vibrating grid . We 
have here focused on the missing link in this chain, the 
slow surface-mediated precursor to rapid turbulence. If 
no other mechanism of vortex multiplication intervenes, 
as is the case in steady flow of 3 He-B with clean smooth 
surfaces, then this process will take over. 



respect to the flow |4|, but it starts from reconnection 
kinks |5|. It is these sharp kinks which prove to be es- 
sential in the simulations for continued generation of new 
vortices . The kinks are primarily produced when an 
expanding Kelvin wave hits the wall. In a rotating cylin- 
der this early stage is marked by roughly linear growth 
in N, as seen in Fig. [3] Its duration in time is similar 
to the measurements in Fig. [2] and it also ends in rapid 
turbulence. The model of this single-vortex instability 
at the wall is studied in Fig. Here we see the expan- 
sion of Kelvin waves from a reconnection kink at the wall 
and a later reconnection again with the wall, which pro- 
duces one new loop. The probability of this process has 
been calculated in the insert. It displays a similar abrupt 
temperature dependence as measured in Fig. ^ 

In numerical calculations self-sustained growth of N 
is not started as readily as in our experiment. Uniform 
flow along a plane wall in Fig. 0] does not support con- 
tinuous growth. In rotating flow the initial configura- 
tion in Fig. |3 had to be specially engineered since exper- 
imentally relevant configurations do not necessarily start 
a continuous process in simulations. Clearly the expla- 
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